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Dispersionless bands, such as Landau levels, serve as a good starting point for obtaining interesting correlated 
states when interactions are added. With this motivation in mind, we study a variety of dispersionless ("flat") 
band structures that arise in tight-binding Hamiltonians defined on hexagonal and kagome lattices with staggered 
fluxes. The flat bands and their neighboring dispersing bands have several notable features: (a) flat bands can 
be isolated from other bands by breaking time-reversal symmetry, allowing for an extensive degeneracy when 
these bands are partially filled; (b) an isolated flat band corresponds to a critical point between regimes where the 
band is electron-like or hole-like, with an anomalous Hall conductance that changes sign across the transition; 
(c) when the gap between a flat band and two neighboring bands closes, the system is described by a single spin- 
1 conical-like spectrum, extending to higher angular momentum the spin- 1/2 Dirac-like spectra in topological 
insulators and graphene; (d) some configurations of parameters admit two isolated parallel flat bands, raising the 
possibility of exotic "heavy excitons" and (e) we find that the Chern number of the flat bands, in all instances 
that we study here, is zero. 



I. INTRODUCTION 

One of the reasons why dispersionless (or flat) bands are 
interesting is that they accommodate, when partially filled, an 
exponentially large number of states. This macroscopic de- 
generacy can be lifted when interactions are added, often lead- 
ing to rich strongly correlated phenomena. The best known 
example is the fractional quantum Hall effect, which arises 
from the degeneracy within flat Landau bands for particles in 
a magnetic field. 

In addition to the Landau problem, other models with flat 
bands have been studied at least since the 1970s, such as amor- 
phous semiconductors^"— . This system is idealized by a lat- 
tice made up of clusters of sites. Both inter-cluster and intra- 
cluster hoppings are allowed. From a mathematical point of 
view, it turns out that the intra-cluster hopping term in the 
Hamiltonian is a projection operator, leading to the existence 
of flat bands^. 

In the 1980s, flat bands were studied in relation to the 
Nielson-Ninomiya theorem^ by Dagotto et alL. They showed 
that it is possible to escape the fermion-doubling problem at 
the price of having an extra flat band in the spectrum. In this 
way, the low energy degrees of freedom of the theory can be 
described as a single Weyl species. 

More recently, Ohgushi et a& studied flux phases in the 
kagome lattice, which is the planar section of ferromagnetic 
textured pyrochlores. If the flux is staggered then electrons 
accumulate a spin Berry phase as they hop. This system con- 
tains isolated flat bands, i.e., they are protected by a gap. On 
the other hand, Bergman et al have studied flat bands with- 
out a gap in similar lattices^ not threaded by fluxes. In their 
models, a flat band is degenerate with one or more other bands 
at a single point, and the touching is topologically protected. 
Each of the works above has identified interesting, but seem- 
ingly disconnected, properties of flat bands. 

The purpose of this paper is to understand the different 
types of flat band spectra, the conditions to obtain them, and 
the properties that follow. This paper is organized around five 




FIG. 1 : Energy dispersions with staggered flux phases <j}+ and 
on the up and down triangles of the kagome lattice . The dispersion 
on the left (type I) corresponds to — 2n and 0_ = — vr while the 
dispersion on the right (type II) corresponds to = = 3n/2. 



main findings. First, flat bands can be isolated by breaking 
time-reversal symmetry (TRS). Second, isolated flat bands 
can be viewed as critical points. On either side of the 
phase boundary the flat band becomes positively or negatively 
curved, which corresponds to transitioning from a particle- 
like to a hole-like band. We also find an anomalous Hall ef- 
fect on either side of the transition, whose sign depends on 
whether the band is electron-like or hole-like. Third, when 
the gap between a flat band and its neighboring bands closes 
with the flat band in the middle, the system is described by 
a single spin-1 conical-like spectrum. This extends to higher 
angular momentum the spin- 1/2 Dirac-like spectra in topolog- 
ical insulators and in graphene. Fourth, one can obtain mul- 
tiple (parallel) flat bands and we provide a concrete example 
of this case. Fifth, we find, for all examples studied here, that 
the Chern number of the flat bands is zero. So as opposed to 
Landau levels, we get no quantized Hall conductance for these 
particular flat bands. 

The model with which we will start our analysis is a simple 
kagome lattice with a tight-binding interaction and two stag- 
gered fluxes, and on alternating triangles. Depending 
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on the values of (j)±, dispersions can be classified into three 
types; (I) a flat band touches two linearly dispersing bands at 
the same point, where the linear bands are reminiscent of a 
"Dirac-like" point, but with spin-1 behavior, (II) an isolated 
flat band that is separated from bands above and below by a 
gap, and (III) a gapless flat band that touches a single massive 
energy band either above or below. The energy dispersions 
corresponding to types I and II are plotted in Fig. [T] Type III 
has been discussed recently^— in hexagonal and kagome lat- 
tices without magnetic flux, and it was found that in this case 
the zero gap is protected by topological arguments'*. Type III 
has been also shown to exhibit a topological insulator phase 
in the presence of spin-orbit interactions^^. Type II appears 
when electrons accumulate a spin Berry phase as they hop 
from site to site, which is equivalent to both up and down tri- 
angles with the same magnetic flus^. Type I necessitates the 
staggered fluxes 0+ 7^ 0_ that we analyze below. 

We will show that the condition for a flat band to occur 
at E — is 0+ + 0„ = TT (mod 2tt). By changing the 
value of the fluxes in such a way that their sum differs slightly 
from TT, the flat band acquires a small curvature, which can be 
positive or negative depending on the values of . Interest- 
ingly, the band curvature implies that if the Fermi energy is 
chosen to be zero, then by tuning the fluxes it is possible to 
change the center of the band from an electron-like pocket to 
a hole-like pocket. This leads to an inversion of the sign of 
the anomalous Hall responsai^. Therefore, the flat band con- 
dition (j)-^ + (j>^ = TT (mod 27r) represents a quantum critical 
point separating two regions with different anomalous Hall re- 
sponses. 

Type I is remarkable in that it displays linearly dispersing 
modes, akin to those of graphene, but differing in two impor- 
tant ways. First, the conical points do not appear in pairs as in 
graphene, but instead there is only one such point within the 
first Brillouin zone (BZ) (see Fig. [T]). Second, these are not 
Dirac fermions (this is why it is possible to evade the doubling 
problem), but instead the effective Hamiltonian in momentum 
space is of the form H = vp k ■ L, where L is the spin- 1 
angular momentum operator (vp is the Fermi velocity). The 
spin-1 -type spectrum (like the spin- 1/2 Dirac-type spectrum 
of graphene) can be viewed as a single quantum spin in a mag- 
netic field, as in Berry's original work on quantum phasesii, 
but with the wave vector k playing the role of the magnetic 
field. Type I does not require the breaking of TRS. However, 
we show that a gap can be opened while leaving the flat band 
untouched by breaking TRS, and thus type I is continuously 
connected to type II . In type II the degenerate states within 
the flat band are protected by the gap at finite temperature, 
and would provide a fertile base to construct correlated states. 

Finally, we will generalize the above results on the kagome 
lattice to other lattices by a formulation similar to Straleyi, 
but again adding staggered fluxes. As an example, we take 
a hexagonal network, similar to graphene, but with three de- 
grees of freedom at each site. We will consider two hopping 
strengths in this case, one for inter-site hopping between near- 
est neighbors (t) and one on each vertex for intra-site permu- 
tations between the three species (g). Each permutation will 
be associated with a flux (f>± on the two sublattices of the hon- 




FIG. 2: The kagome lattice with fluxes (j>+ and (j)- on alternating tri- 
angles, which correspond to a flux —((/>+ + (/)_) inside each hexagon. 
si,2,3 are vectors pointing from the center of a down triangle to its 
three up neighbors. 



eycomb lattice. We will show that the kagome model can be 
obtained from this in the limit t ^ g. Additionally, we will 
show that the honeycomb lattice admits two parallel flat bands 
that are isolated from each other and all other bands, which is 
a feature that the kagome lattice does not have. 



II. FLAT ZERO-MODE BAND IN THE STAGGERED-FLUX 
KAGOME LATTICE 

Consider the tight-binding Hamiltonian defined on a 
kagome lattice, where staggered fluxes cf)^ and (/)_ are ap- 
plied within alternating triangles ("up" and "down" triangles, 
respectively), as shown in Fig. [J] 

For convenience define the phase factors a± = e*"*'±/3. Let 
si = (0,-1), S2 = (V3/2,l/2), and S3 = (-^3/2, 1/2) 
be the vectors pointing from the centers of an up triangle to 
its three down neighbors, and define ~ e-*k (si-sj )^ ssixHa 
j = 1,2,3. In momentum space, the Hamiltonian can be 
written as: 

(0 a+ + a^dl^ + ct-d]^ \ 

+ + ,(1) 

a+ + a^dl^ g / 

where g is the hopping strength. The characteristic polyno- 



3 



mial for this matrix is 

PiE)^-E''+g^ai(k)E + g^ao{-k) , (2) 

where 

ai(k) = [3 + a+a_ g(k)] + c.c. 

ao(k) = [(a+ + (a^Q- +a+a_^)g(k)] +C.C. , 



and g(k) 



"k 



j23 
"k 



The condition for a flat (k- 



independent) band is obtained by setting the overall factor of 
q(k) in P{E) to zero. 



which is equivalent to 

E = -g{a\ + S-^) . 







(3) 



(4) 



Combining Eqs. (|2]i-(|4|i gives the following equation for the 
energy eigenvalues of the flat bands: 



E{E - V) = , 
which has three possible solutions: 



(5) 



(a) E ^ flat band: this configuration is achieved for 

(/)+ + 0- = TT (mod 27r), 

(b) E = —2g flat band: this configurations is achieved for 
= 2Tm±, with n± integer valued numbers, and 

(c) E = 2g flat band: this configurations is achieved for 
'P± = 7r(27i± + 1), with n± integer valued numbers. 

The cases E = ±25 (type III) are similar to those discussed 
in RefS]^— , where the flat band touches a parabolic electron- 
like band at its bottom (for E = —2\g\) or a hole-like band at 
its top (for E = +2\g\). Here, we shall focus instead in the 
case where the flat band is at = 0. In Fig. [T] (types I and 
II) we show two particular choices for (p±, which are repre- 
sentative of what we classify as types I and II spectra. Type 
I contains a cone vertex touching at the point k = 0, which 
we illustrate by setting (f)-^- = 3n and (this choice of 

phase can be interpreted as tight-binding hoppings —g for up 
triangles and +g for down triangles). Type II contains an iso- 
lated flat band, which we illustrate by setting <f>± = 3tt/2 
(this choice can be interpreted as tight-binding matrix ele- 
ments ztig for hopping anti-clockwise or clock- wise around 
the triangles). Notice that because triJk = and one of the 
eigenvalues is E = 0, the other two eigenenergies must sat- 
isfy E+ (k) + (k) = 0, so the spectrum is symmetric with 
respect to zero in both types I and II. 

The general condition for nodal touching (type I) can be ob- 
tained by requiring that there is another E = eigenvalue, so 
that at least one other band touches the flat band. When such a 
solution exists, the derivative of the characteristic polynomial 
P'{E) also has a zero at E' = for some value of k. This 
condition translates to 



which admits three different solutions: 

(A) Nodal point at F = (0, 0): this type I configuration is 
obtained if the condition 0+ — (/)_ = Stt + Q-kti is satisfied, 
where n is an integer, and it is illustrated in Fig. [T](type I), 



(B) Nodal point at k = 
uration is obtained if the condition 
satisfied, where n is an integer, and 



3v^ 

+ 



, 0): this type I config- 
— (b- = Stt + 67rn is 



(C) Nodal point at k = = (^'S^'O)- ^^P^ ^ 
configuration is obtained if the condition (f)+ — (p- = tt + 67rn 
is satisfied, where n is an integer 

Even though types I and II are particle-hole symmetric with 
a flat band at E' = and have staggered fluxes obeying the 
constraint (/)++(/)_ = tt (mod 2tt), type II lacks the nodal 
conditions (A)-(C) mentioned above. 



A. Nodal touching and the spin-1 cone 

Let us now expand the Hamiltonian, in type I , near the 
vertex point for small |k| = 



At the same time 



we move into type II by applying a slight flux offset from the 
condition for the touching: 0+ = 3(7r + 5) and 0^ = —3(5. 
We will interpret 5 as a "mass" term. To first order in 5 and k 
the Hamiltonian becomes: 



L'^ + ky L'y + 2\I-SL: 



= g-j= {kx,ky,m) ■ L' 



(7) 



where m = 2-^/2/3 d. It is straightforward to check that the 
matrices 



(8a) 




ai(k) = ^ a+a_(7(k) 



(6) 



(8b) 



(8c) 



satisfy the angular momenta algebra [L'^,Ly] = iL'^ (along 
with the cyclic permutations of x, y, and z), and that they have 
eigenvalues —1, 0, +1, i.e., they form a spin-1 representation 
of SU(2). 

The eigenvalues of the Hamiltonian (|7]i are i?k = 

g-^ Jkl + kl+m? where £k = -1,0, +1 is 
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FIG. 3: Left (right): middle band showing electron-like (hole-like) 
dispersion corresponding to e < 0(> 0) in Eq. ilOl . 

the eigenvalue of angular momentum along the direction 
{kx, ky, m). Therefore we obtain the three bands, with the flat 
band being the one with zero angular momentum. The other 
two bands describe the cone when (5 = 0, and two parabolic 
bands separated from the flat band by a gap A = g 2a/3 5 
when 5 is non-zero. (A particular instance of the 5 = spec- 
trum has recently also been predicted in Ts optical lattices by 
Bercioux et al.—) 

The spin-1 structure has interesting topological properties, 
namely it can be viewed as a generalization of the Berry phase 
for the spin- 1/2 spectrum of graphene. One implication is that, 
when the gap is open by breaking TRS, the upper and lower 
bands have a quantized Hall conductance. In the Appendix we 
explicitly compute the Chern number over the first Brillouin 
zone for the three bands. 



B. Transitioning between electron and hole bands 

We now consider deviations from the flat band condition 
for the case when the middle band is isolated, as in type II . 
For concreteness, consider the case </>+ =0- = 3(7r/2 — e). 
For small e, the energy of the middle band will be close to 
i? = 0, so we can obtain the dispersion for the middle band 
by dropping the cubic term in the characteristic polynomial 

E{\.) = -g^ + 0{e'). (9) 
ai(k) 

Expanding ao(k), ai (k) up to order |kp and e, one obtains 

Eik)=g (^4e-||kp^ , (10) 

corresponding to a band mass 771^ = — 3e/2 g for the middle 
band, so it has a hole-like dispersion for e > and an electron- 
like dispersion for e < 0, as depicted in Fig. [3] This trivial 
mathematical result is physically remarkable in that one can, 
in principle, change the character of a band from electron-like 
to hole-like by varying one parameter in the Hamiltonian. 

Notice that there is also a band shift 4eg, which adds to 
the chemical potential. There is an interesting result when the 
chemical potential is fixed to fi — 0: the Fermi surface is 



pinned and independent of e. It is best to see this effect prior 
to any perturbation in e or expansion in | k| . The Fermi surface 
is in this case the locus of k points for which P{E = 0) = 0, 
those that satisfy ao(k) = 0. For the electron-like case, the 
Fermi sea is in the region bounded by the ao(k) = sur- 
face that contains the F-point, whereas for the hole-like case 
the Fermi sea is the complementary region in the Brillouin 
zone. The anomalous Hall effect is given by the integral of 
the Berry curvature over the Fermi sea. As we show in the 
Appendix, the Chern number, the Berry curvature integrated 
over the complete Brillouin zone, for the middle flat band is 
zero. This means that the sum of the anomalous Hall effects 
for the electron-like and hole-like Fermi seas is zero. Thus, 
as one tunes across holding /i = 0, the anomalous Hall effect 
will change sign. Of course, by tuning fi one can vary the 
anomalous Hall effect continuously. 

C. Time-reversal and particle-hole symmetries 

First, we discuss time-reversal symmetry. Consider a tight- 
binding model of spinless fermions described by 

H = ^ V4(k)i?.,(k)^j (k) , (11) 

k 

where ip is an annihilation fertnionic operator and k takes val- 
ues on the first Brillouin zone. Under time-reversal transfor- 
mation, 

^ ^ ^V4(-k)JI*(k)^,(-k) 

k 

= ^V4(k)i^* (-k)^,(k) . (12) 

k 

For Ti to be time-reversal invariant, one way would be to have 
i/(k) = k). Looking at our Hamiltonian more care- 

fully, though, we see that there is a freedom to redefine the 
hopping matrix elements without changing the fluxes <j>^ and 
This freedom can be mathematically described by the 
following gauge transformation: 

H{k) H{k) = AH{k)A^ , (13) 

where 

/e*"i \ 
A= e^"2 . (14) 
\ e'"3 J 

Hamiltonians H{k) and H{k) related by the gauge transfor- 
mation Eq. ( fTSl ) represent physically equivalent descriptions 
of the system. The requirement of time-reversal symmetry, 
taking into account the gauge invariance given by Eq. ( fT3] ), 
becomes then 

H{k) ^ AH*{-k)A'< , (15) 
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from which we get the following conditions 



^{2i/3)4>± _ 1(01-02) _ „i(a2-a3) _ 1(03-0:1) 



(16) 



Equation ( fTSI l immediately implies that symmetry under time 
reversal is satisfied if 0+ = n+7r and (f>^ = n^n, for integers 
n+ and n_ such that n+ — n_ = 3/ (I integer). This is exactly 
equivalent to the condition for (t>± such that the spectrum has 
a gapless flat band at ii^ = 0, i.e., type 1 . Therefore, we con- 
clude, for this given model, that in order to have an isolated 
flat band time-reversal symmetry must be broken. 

Now consider particle-hole symmetry. The key observation 
is the following: when the spectrum of H has three bands, 
as in the kagome lattice, particle-hole symmetry only exists 
when there is a flat band at = and the two other bands 
have opposite energies. We have already worked out the con- 
ditions for the existence of a flat E = Q band in the kagome 



lattice to be 



+ 0- 



TT, which also dictates the conditions 



for particle-hole symmetry (if the Hamiltonian is a traceless 
matrix). Notice that T-L can be particle-hole symmetric with- 
out being time-reversal invariant. (A spin-1 cone is a situation 
where both symmetries are present.) We will use this impor- 
tant aspect later when we calculate the Chern number of the 
bands. 



III. HEXAGONAL LATTICE MODEL 

Having discussed the main properties of the kagome model, 
we now consider an alternative model defined on the hexago- 
nal lattice. We will show that these models are closely related 
and that the kagome model is a limiting case of the hexagonal 
one. Even more interestingly the spectrum of the hexagonal 
model contains two parallel flat bands that are separated from 
each other and from all other bands by a gap. 



A. Deflnition of the model 

Consider a tight-binding model on a hexagonal lattice, 
where the particles have three flavors. Define the six- 
dimensional basis of particle operators ^), where 
^jl — 1, 2, 3 is the flavor index, and a, b are the two sublattices. 
In real space the Hamiltonian is: 

^ = ^ E ^Ia. '^Mli-b - r„) i'b^,, + H.c. (17) 

{ah) 

+.9 X] '^i,ti.BtJ.ua+ipa,u + ipl^^Bfj,^a-i:b,u + H.c. 

The first term (coupling t) is a nearest-neighbor hopping be- 
tween the two sublattices that conserves the flavor index, but 
has correlated flavor-direction hopping controlled by h^. 



/lM(r(, -Va)=\ 

0, otherwise. 



(18) 



The second term (coupling g) is the on-site phase dynamics, 
where B is the 3x3 cyclical permutation matrix. 



B 




B' =B^ 



B-' 



(19) 



and the on-site flavor changing phase factors are a± = e"'^±/'^ 
(using the notation of Sec. II). In other words, inter-site hop- 
ping conserves the flavor index. Flavors can be permuted on- 
site, and each permutation is accompanied by an on-site phase 
factor 0±. This Hamiltonian can be diagonalized analytically 
for particular values of (j)± . 



B. Equivalence to the kagome model 

Let us write the wave function basis as 

W'l,i,V'M''/'l2'V'l2'V'I,3'V'l,3) and work in momen- 
tum space. In this basis, our Hamiltonian takes the following 
six-dimensional form: 



Ti G 

T2 G 
G Gt Ta 



(20) 



where 







G 



(21) 



In the extreme case when g ~ Q, H \s block diagonal in the 
T matrices and is trivial. It has two triply degenerate eigenval- 
ues E = ±t and eigenvectors vj^ _|_(k) = (1, ±e^^'^^') / \/2. 
This corresponds to a "dimerized" state since the particles 
cannot hop between sublattices. (Recall that hoping between 
sublattices is directionally controlled by the flavor, so if the 
particles are not allowed to change flavor in the vertices be- 
cause g = 0, they can only hop back and forth within a given 
bond between an a and a b site.) 

Now project the full Hamiltonian to the lowest energy 
states, E = —t, and expand to first order in g / 1. The projected 
Hamiltonian, H, is given by the matrix elements if^^ = 



For example, _ffi2(k) = + dl^a^)/\/2. 



and so on. The result is that i/(k) is identical to the Hy^ on 
the kagome lattice in Sec. II, up to an overall additive energy 
-t. 

To check the equivalence explicitly, let us solve the Hamil- 
tonian at one particular choice of fluxes: (/)+ = (/)_= 0, for 
example. In this case the spectrum is given by 



E 
E 



±t 
9 



9 



(22) 



1 

2±2 



4^2 + 9(?2 ± 4t5V3 + g(k) + g(k) .(23) 



6 




FIG. 4: Exact spectrum for t — 1, g = V3/2 and 4>± = 3tv/2. 

To first order in g/t the eigenvalues above Eq. ( |23] | reduce to 

£ = -t-.g (24) 
S = -i + |(l ± v/3 + g(k) + g(k)] , (25) 

and similarly with t —t. This is exactly the spectrum of 
a kagome tight-binding model in zero field, or type III in our 
nomenclature, confirming that the kagome model is a limit of 
the hexagonal model. Note that in this particular case, where 
time-reversal symmetry is not broken, the flat bands touch 
parabolically at k = 0, as discussed before in Sec. II. In fact, 
in the exact solution in Eq. (l23T l for all values of g and t the 
flat bands are not gapped for the same reasor>i^. 



C. Parallel flat bands 

There is a rather interesting configuration of parameters 
for which the spectrum contains two distinct flat bands. As 
in the kagome model, breaking of time-reversal symmetry is 
required to create gapped flat bands. We thus break time- 
reversal symmetry and choose the phases (/>+=(/>_= 37r/2. 
Further, we tune the interactions such that g = (V3/2) t. The 
spectrum for this Hamiltonian is shown in Fig. |4] Notice 
that, remarkably, because of the two flat bands even non-zero 
momenta particle-hole excitations in this system would have 
macroscopic degeneracy, and exotic heavy excitons could be 
formed. 

We close this section by noting that there is at least one 
other hexagonal model where this parallel band structure can 
be induced by breaking time-reversal symmetry. Wu and 
Das Sarma recently studied ground-state properties of inter- 
acting spinless fermions in the p^.j,-orbital bands in the two- 
dimensional honeycomb optical latticei^. They considered a 
CT-bonding interaction, which describes hopping between p or- 
bitals on neighboring sites when the orbitals are oriented along 
the bond direction. They found gapless flat bands since their 
model does not break time-reversal symmetry. However, we 



would like to point out that one can add to their model a time- 
reversal symmetry breaking interaction, on site, between the 
Px and Py orbitals of the form gi {p\py — vXvx)- Just as in 
our example, parallel flat bands appear when g is tuned to a 
special multiple of the hopping strength t. 

IV. SUMMARY 

Dispersionless bands can be the starting point for construct- 
ing strongly correlated electronic states. The lack of electron 
kinetic energy leads to a macroscopic degeneracy when dis- 
persionless bands are partially filled, and interactions become 
responsible for lifting the degeneracy and selecting the many- 
body ground state. This situation is the case for Landau levels, 
which are flat bands created by an external magnetic field. 

In this paper, we have analyzed different types of spectra 
that contain flat bands in tight-binding systems in the pres- 
ence of staggered fluxes. We have seen that it is possible to 
separate a flat band from the other bands by a gap when time- 
reversal symmetry is broken. In these situations, the flat band 
can be viewed as a critical point (with zero curvature) that sep- 
arates electron-like from hole-like bands, and we can switch 
between these two curvatures by changing parameters in the 
Hamiltonian. 

When the gap is closed and the flat band lies between two 
other bands, one obtains a spin-1 conical spectrum, extending 
to higher angular momentum the spin- 1/2 Dirac-like spectra 
in topological insulators and in graphene. 

We have also presented examples of tight-binding systems 
where it is possible to obtain more than one isolated flat band. 
Specifically, we showed examples with two isolated parallel 
flat bands. 

Although we made progress in understanding several as- 
pects of flat bands, two points remain open questions and de- 
serve further investigation. First, we do not have a generic 
proof that time-reversal symmetry must be broken to isolate 
a flat band. Nonetheless, it is natural to speculate that this is 
true in general, as it holds in all examples that we have found, 
in addition to the well-known case of Landau levels. And, 
second, in all examples discussed here as well as in a class of 
tight-binding models defined on a line graphic, the flat bands 
have zero Chern number Whether this is an intrinsic property 
of these bands remains unclear to us. If it is possible, however, 
to find examples of flat bands with non-zero Chern number in 
the absence of an external magnetic field, this could be an in- 
teresting scenario for realizing strongly correlated electronic 
states with topological order. 

Note added in proof. Recently, we learned of a promising 
realization of flat band systems by Koch et al. using circuit- 
QED based photon latticesi^ 
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Appendix A: Chern numbers for bands in the staggered flux 
system 

Haldane^ showed in a seminal paper that it is possible for a 
system to exhibit the quantum Hall effect without Landau lev- 
els provided the system breaks TRS. Recently, in the context 
of the anomalous quantum Hall effect, Ohgushi et afi pro- 
posed a three-band model with electrons hopping in a kagome 
lattice in the presence of a background spin texture. In their 
model, the spin texture opens a gap in the spectrum and gives 
rise to a Berry phase such that the Chern numbers of the bands 
are 1, 0, —1 when TRS is broken. In this appendix we com- 
pute explicitly the Chern numbers in the type II spectrum of 
our kagome model. 

Recall that when (/)++(/)_ = tt and 0+ — 0_ = Stt we 
have a flat band with a Dirac point as shown in Fig. [T] For 
this choice of fluxes, as discussed previously, the Hamiltonian 
is time-reversal invariant and the system presents no Hall re- 
sponse (type I). However, when the flat band is maintained but 
a gap is opened, time-reversal invariance is lost and we have 
the possibility of bands with non-zero Chern numbers (type 
II). We parametrize the gap by a mass term 6, such that the 
fluxes 0+ = 2tt + 6 and (/)_ = — tt — (5. 5 = corresponds 
to a Dirac cone at the center of the BZ. For reference we 
give the complete energy spectrum, although we will expand 
around small k below: E± = ±-\/ /(k) and Eq = 0, where 
/(k) = 6 — 2 j cos [k • {sj — Si) — 26/3], and the sum- 
mation is over the cyclic permutations {i,j) = (1, 2), (2, 3), 
and (3, 2). 

The Chern number of the nth band is defined as the sum- 
mation over the first BZ, 



{nk\Xjc\mk){mk\Jy\nk) - (J^; ^ Jy) 



iEnik)-EmiW^ 



^ Vk X A„(k) 

keBZ 



(Al) 



keBZ 



Here, i3„(k) is the field strength associated with the Berry 
vector field A„(k) = — i(nk| Vk|n.k) and J = {Jx,Jy) is 
the current operator given by J = W\^H. One can see that 
J2n^n = by the antisymmetry of C„ as x and y are inter- 
changed. Because the Chern numbers of the bands are topo- 



logical quantities-^^^, their values can only change when a 
band touching occurs. Choosing 6 to be very small, an ar- 
bitrarily small gap to ~ 5 is opened (still keeping the flat 
band) and a near degeneracy appears for k w 0. Around 
this point, the Hamiltonian is that of a spin-1 system with 
iJk ~ kxLx + kyLy + mLr. 

To perform the summation explicitly it is conve- 
nient to define the vector / = {kx,ky,m) = 
|/|(sin0 cos (/), sin sin0, cos 0), which can be viewed as a 
magnetic field in momentum space coupled to the spin oper- 
ator (c.f.. Berry''*). First, we compute the eigenvectors of iJk 
with the respect to the z axis and then we apply a rotation 
to bring the spin states to an arbitrary {9, 0) direction. Let 
{|x+)7 Ixo), Ix-)} be the eigenstates of with eigenvalues 
1, 0, —1 respectively. The eigenvectors of _ffk in a general di- 
rection (9, (p) are given by 



\lpn 



-iOL-u 



(A2) 



with n = +, — , 0, which in explicit form reads 



( 1+cos \ sin 6 ^icp f 1 — cos 6 \ \ / A a \ 



IV'o) 



= I 



and -ffk| 0: 



V2 



^ — i(t> sin 6 rT\c n ^i(p sinO 

^ , cos p, e ^ 



(A3b) 



-i<f> 



sin 9 
' V2 ' 



(1+^) ) ^ (A3c) 



±^|k|2+TO2|V^±) and Hi^lijo) 



0. A 



straightforward calculation of the field strengths gives us 



B±(k) = ± 



(to2 + |k|2)3/2 



Bo{k) = 



(A4) 



The contributions of these fluxes to the Chern numbers are 
found to be ±sgn(7Ti) and zero. We have also confirmed 
this result numerically over the entire BZ without linearizing 
around the F point. As we cross the gap, 5 (equivalently, to) 
changes sign and the Chern numbers of the upper and lower 
bands change sign as well, while the Chern number of the flat 
band remains zero. Because the topologolical nature of the 
Chern numbers, their values will remain unaltered until a new 
band touching occurs, which will happen for 5 = ±n, when 
the Dirac point moves to one of the corners of the BZ. 
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